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ABSTRACT 

We apply the hierarchical clustering model and non-linear perturbation theory to 
the cosmological density and temperature fields. This allows us to calculate the in- 
tergalactic gas pressure power spectrum, SZ anisotropy power spectrum, skewness and 
related statistics. Then we show the effect of the non-gravitational heating. Our model 
confirms recent simulations yielding mass weighted gas temperature Tg ~ 0.35 keV and 
reproduces the power spectra found in these simulations. 

While the SZ effect contains only angular information, we show that it is possible to 
extract the full time resolved gas pressure power spectrum when combined with galaxy 
photometric redshift surveys by using a variation on the cross-correlation. This method 
further allows the disentanglement of the gravitational and the non-gravitational heat- 
ing. 

Subject headings: Large scale structure-cluster-cosmology; CMB; Perturbation theory 



1. Introduction 

One of the open questions in cosmology today is the state of the intergalactic medium (IGM). 
In the present popular cosmological models, a putative cosmological constant accounts for perhaps 
2/3 of the energy density of the universe, and a similarly mysterious cold dark matter component 
accounts for 2/3 of the remainder. Of order 10% of the energy density is supposed to be accounted 
for by baryons (Lange et al. 2000; Tytler et al. 2000) as obtained from primordial nucleosynthesis 
calculation. Surprisingly, the vast majority of the baryons in the local universe is as of yet unde- 
tected. The directly observed known components consisting of stars, as well as hot and cold gas in 
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galaxies, groups and clusters, only account for a few percent of this baryon budget (Fukugita et al. 
1997). Cluster of galaxies are known to possess a large gas fraction, of order 20%/i5q (Pen 1997). 
This gas fraction, if representative of the universe at large, is sufficient to account for the majority 
of baryons in the universe (White et al. 1993). Since gravity acts on all particles equally, one 
expects this ratio within the virial radius of clusters to be representative of the universal fraction. 
In the standard bottom-up picture of structure formation, clusters are the most massive objects 
in the universe, and therefore the latest to form. This suggests that the majority of baryons is 
in a diffuse intergalactic form, which forms the intracluster medium when clusters collapse. If the 
baryons were in other forms, for example compact objects which formed early in the universe, it 
would be extremely difficult to understand how those baryons would be released back into diffuse 
form when clusters form. 

Pen (1999) showed that the IGM must have been heated by sources other than the gravitational 
energy of collapse and its resulting shock waves. Its direct detection would be of significant interest 
to quantitatively understand the distribution and state of the IGM, and thereby infer its thermal 
history. Unfortunately, hot diffuse gas is difficult to detect in emission. Luckily, the Cosmic 
Microwave Background (CMB) scatters off all free electrons through inverse compton scattering, 
allowing us to "see" the IGM through scattering, in analogy to absorption. This is known as the 
Sunyaev-Zel'dovich effect (Zel'dovich and Sunyaev 1969), which is becoming routinely observable. 
It becomes the dominant source of CMB fluctuations on arc-minute scales. Several experiments to 
conduct deep blank field or all sky search are under development. The upcoming CMB experiments 
such as AMIBA (Array for Microwave Background Anisotropy (2003)), South Pole Submillimeter 
Telescope (2003), Planck (2007), etc. and the current CBI (Cosmic Background Imager (2000)) are 
capable of detecting this angular scale. Using multi-frequency information, the SZ fiuctuation can 
be disentangled from the primary anisotropics at all angular scales (Cooray et al. 2000a). 

On the theory part, various approaches have been explored to compute the SZ angular power 
spectrum. Amongst these, the Press-Schechter formalism is perhaps the most widely used due to 
its versatility and ease of implementation (Cole and Kaiser 1988; Makino and Suto 1993; Artio- 
Barandela and Mucket 1999; Komatsu and Kitayama 1999; Cooray 2000b). The authors above 
adopted the cluster gas model, while Da Silva et al. (1999); Refregier et al. (1999); Seljak et 
al. (2000) used simulations and Cooray et al. (2000a) used a simplified gas pressure bias model 
motivated from simulations to probe the statistics of the SZ effect. 

One of the shortcomings of the SZ effect is its lack of redshift information. Since scattering is 
independent of redshift, the SZ effect on one hand allows a direct probe of the IGM to high redshift, 
on the other hand makes it challenging to disentangle the contributions arising from different 
redshifts. In this paper, we apply the hierarchical clustering model and non-linear perturbation 
theory to directly compute the SZ effect. This analytical method enables us to to extract distance 
information of intergalactic gas from the cross correlation of the SZ effect with galaxy surveys and 
at the same time, check the results of both simulations and the Press-Schechter formalism. In 
section 2, we develop our gas pressure model and in section 3 the SZ statistics including power 



-3- 



spectrum and bispectrum. Section 4 contains our results on the SZ-galaxy cross correlation and 
the method to extract the redshift distribution of the SZ effect. The effect of non-gravitational 
heating and the method to extract it from overall gas pressure power spectrum are discussed in 
section 5. We discuss the potential inaccuracies arising from the approximations in section 6. The 
paper concludes with section 7. 



2. Gas pressure power spectrum 

The temperature distortion caused by the SZ effect (Zel'dovich and Sunyaev 1969) is: 

e(n) ^ ^^ill^ = [4 - x/tanh(x/2)] ^ -2yS{x) (1) 

J^CMB - 1 

where n is the direction on the sky, x = hu / (kTcMB)-> and the scattering function S{x) — > 1 when 
x <C 1 (Rayleigh- Jeans tail). In this limit the SZ effect results in an apparent cooling of the CMB 
background. The "y" parameter is defined as 

y{n) = / UekTgdl = — 4 / Pe{n)dl = / Peyp{n)dl (2) 

Here, Tg and Ue arc the temperature and number density of free electrons, respectively. Pg is the 
gas pressure, yp = Pe/ Pe = + ^g)Tg/Tg and Tg = {{1 + 5 g)Tg) is the gas density weighted mean 
temperature, dl = a{z)C{x)dx is the proper distance, a{z) is the scale factor, x{z) is the comoving 
distance. 



C{x) = ^ (3) 

describes the geometric effect of the curved universe, K = —1, 0, 1 for open, flat and closed universes 
respectively, and Rq = ■^(l — Qo)~^^^ is the curvature radius, po = Pc^o is the present cosmological 
matter density. 

Pen (1999) showed that the intergalactic medium has most likely been preheated by non- 
gravitational sources to ~ 1 keV per nucleon in order to be consistent with the observed upper 
bounds from the X-ray background. We adopt the model of (Pen 1999) to express the gas distri- 
bution as a convolution of the matter distribution: 

5g{x) = j 5{x')Wg{\x-x'\)d^x' (4) 

This equation expresses gas as being less clumped than dark matter, partly due to the required 
preheating. The effective radius in the top-hat window function, which is the gas heating radius, 
has the typical value ~ lh~^ Mpc from the X-ray background constraint (Pen 1999). For simplicity, 
we adopt the Gaussian window Wg{r) = exp(— r^/2rg)/(\/27rrg)^ {vg ~ 1/3/i^^Mpc corresponds 
to Ihr^ Mpc top hat window). The window function for specific heating models can be obtained 
from hydrodynamic simulation by comparing the gas power spectrum to the dark matter power 
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spectrum (Ma and Pen, 2000). We also need to relate the gas temperature with the density. First, 
we consider the gravitational heating. We adopt the cosmic energy theorem (Peebles 1980) for 
the gas temperature model, which specifies the ratio between gravitational binding energy and 
total kinetic energy K. The pressure depends on the thermalized fraction of K. The translational 
kinetic energy is thermalized from the energy released when particles shell cross. A model of the 
thermalized energy is thus given by the difference in energy between two particles separated by 
a non-linear scale in Lagrangian space, which is the distance at which they can be expected to 
have shell crossed. The exact procedure amounts to solving the non-linear evolution equations 
directly. But wc can treat the effect statistically in a linear fashion. In the initial linear evolution, 
the gravitational potential remains constant. After virialization, the gravitational energy at a fixed 
location remains almost constant. In an Eulerian description, we can describe the energy of particles 
at a final virialized location as the energy released as a particle travels from its initial position to 
the final virialized location. While the initial position is not exactly known, we take a spherical 
average over the non-linear scale to average over all possible initial locations. We thus have 

* is the gravitational potential, V^* = —AnGpa^d. ^'(x) = / ^{r)We{\x — r\)d^r (hereafter call 
the "electron window function") is the potential averaged over the non-linear scale re- We choose a 
Gaussian window function We(r) with the the non-linear scale rg, which for 2; = is re ~ 5h~^Mpc. 
Hereafter we will adopt this value of rg. X = 0.76 is the mass fraction of the hydrogen in baryonic 
matter. Then, {kTg)k oc 6{k) [I - We{k)] /k^ = 6{k)fe{k). Here, We{k,z) = exp{-k'^rl{z)/2) is the 
Fourier transform of the electron window function. Equation (5) has some unphysical statistical 
properties. The spatial average of the temperature, for example, is exactly zero, and for our 
model using Gaussian random fields, it will be negative in half the volume. But for purposes of 
modeling the SZ effect, temperature is only observable when multiplied by density, and regions of 
positive temperature will have high density, while the negative temperature regions only contribute 
negligibly to the SZ effect. This model is not meant to be an exact description, but hopefully 
captures the statistical properties, while being simple and thus exactly solvable. 

In order to estimate the accuracy of these assumptions, we first compute the gas density 

weighted temperature and the mean y parameter. Wc define S(x) = (27r)~^ J 5^ ex'p{—ik-x)d^k, the 

power spectrum P{k) = and the variance A^(A;) = k^P{k)/2iT^. We adopt the initial power 

spectrum (Peebles 1983; Davis et al. 1985) Punearik) oc A;^+"/(l + 1.7k + 9k^-^ + k'^) (here, k is in 

unit of rio^^/Mpc and we choose the Harrison-Zel'dovich- Peebles scale invariant spectrum n = —1 

corresponding to a = 0), cluster-normalized density fluctuation at 8/i~^Mpc by as = 0.53fi^^ 

[p = 0.53 for ACDM and /? = 0.45 for OCDM ] (Pen 1998) and the Peacock and Dodds (1996) 

fitting formula to convert the linear Punear to the nonlinear P{k). Then 

1 Am u Hk 
f, = {{l + 6,)T,)=AnGpo-j^{l + z)j^ A\k)W,{k)fe{k)- (6) 



rrieC^ 
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Wg{k) is the Fourier transform of the gas window function in equation (4), pB = pc^B is the 
present baryonic matter density and fy is the mass fraction of baryonic matter in gaseous form. 
Since cluster gas fractions (Danos and Pen 1998) are comparable to the baryon fraction obtained 
from big bang nucleosynthesis, we expect that fg ~ 1. We use fg = 0.9, fish? = 0.02 and the 
dimensionless Hubble constant h = 0.67. The results for SCDM, OCDM and ACDM are shown in 
table 1 and figure 1. 

In the absence of non-gravitational heating effects, we find the present Tg is around 0.35 
keV and the mean SZ temperature distortion ~ 5 x 10~^K . They are all consistent with the 
simulations and Press- Schechter formalism results (Refregier et al. 1999; Seljak et al. 2000). 

We define the pressure correlation function ^p(r) = {yp{x)yp{x + r)) = ([1 + Sg{x)]T{x)[l + 
5g{x + r)]T{x + r)) jT^ . ^ Then the pressure power spectrum 



Pp{k,z) = 



I 



A\k,z)Wg{k)fe{k,z) 



dk 
T 



(2vr)f 



54(ki, k2, ka, k4; z)Wg{k^)feik2,z)Wg{ks)fe{k^, z)(fk2d^k4 



(8) 



(27r) 



^ J Wg{ki)P{h,z)fe{k2,z)P{k2,z) [Wg{ki)fe{k2,z)+Wg{k2)feiki,z)]d'k2 



/ B3{ki,k2,k;z)Wg{ki)fe{k2,z)fe{k,z)d^k2 
+P{k,z)f^{k,z)) 



Statistical isotropy implies, ki +k2 = — ka — k4 = k. The first term in the parentheses is from the 
non-Gaussian (connected) part of the four-point correlation {5g{x)T{x)dg{x + r)T{x + r)). B4 is the 
density polyspectrum. The second term is from the Gaussian part of {Sg{x)T{x)Sg{x + r)T{x + r)). 



^Our definition is different to tlie usual definition of ^p(r) = {Sp{x)6p{x + r)). Here, Sp = yp — 1 = (Pe — Pe)/Pe- 
These two definitions only differ by a constant 1 and tlic resulting power spoctrums only differ when A: = by a Dirac 
function. We choose this definition because this is the most convenient way to deal with SZ effect. 



Cosmologies 


SCDM 


OCDM 


ACDM 




1 


0.37 


0.37 




0.53 


0.83 


0.90 


rg/{h-^Mpc) 


1/3 


1/6 


1/3 


1/6 


1/3 


1/6 


fg{z = 0)/keV 


0.3 


0.35 


0.30 


0.35 


0.33 


0.38 


y X 10^ 


1.4 


1.7 


2.9 


3.5 


2.3 


2.6 



Table 1: Present day mass weighted average gas temperature Tg{z = 0), y, and temperature distor- 
tion skewness parameter 63 = ^|^Q_J'^2p for various cosmologies. Two different gas-mass relations 
are used (r^ = 1/3, 1/6), which only has a small effect on the results. 



-6- 



This part also produces a term (27r)^(5/j(k) in Pp{k), which only affects the spatially averaged 
quantities such as (P^) and will be included explicitly only when it is observable. The third term is 
from the three-point correlation {5g{x)T{x)T{x + r)) + permutation. i?3 is the density bispectrum. 
The last term is from the raw temperature correlation {T{x)T{x + r)). 

We adopt the hierarchical model (Fry 1984) which allows all higher order correlations to be 
expressed as the sum of products of two point correlations over all configurations, so 

i?,v(ki,...,kN) =^QiV,a W^AB (9) 

a=l labelings edges 

Pab is the two point power spectrum, a denotes different configurations. Coefficients QN,a generally 
depend on configurations. In the highly non-linear regime, they degenerate (Scoccimarro and 
Prieman 1999). Because the SZ effect is mostly contributed by the highly nonlinear regime, we adopt 

the saturated Qab in the highly non-linear regime: * = ^i+6~xV+3»+i+l"+i ^i"* = 
from HEPT (the Hyper- Extended Perturbation Theory) (Scoccimarro and Prieman 1999). Here, 
n is the linear power spectrum index and we choose the value oi n &i k = '^^i{ki/N). In this 
framework, 

B3{ki,k2,k3) = Q3{PlP2 + P2P3 + P3Pl) 

B4{ki,k2, ks, k4) = Q4 [P1P2P3 + (3 Permutations) -|- P1P12P4 + (H permutations)] (10) 

where Pi = P{ki) and Pij = P(|ki + kj|). 

Results arc shown in figures 2 and 3. Calculations show that the non-Gaussian term is dominant 
in the spectrum even when k is small and this non-Gaussian term itself is mostly contributed 
by the nonlinear regime. Besides the reason that the electron window function sweeps off most 
contribution from the linear regime, these behaviors reflect the strong correlation between S and T 
and the domination of the highly nonlinear (non-Gaussian) regime. This argues that the isotropic 
simplification for Qtv is at least self-consistent. We define the gas pressure bias bp as bp{k, z) = 
Pp{k, z)/P{k, z). The amplitude of bp varies only weakly with time. This is an expected consequence 
of the hierarchical model, argued as followed. The time dependence through the density power 
spectrum is basically canceled because Eq. (8) shows that the dependence of Pp{k, z) over time is 
roughly Pp{z) oc B4{z)/P^{z). The hierarchical model implies that roughly -84(2;) oc P^{z). Then 
bp{z) = Pp{z)/ P{z) oc P^{z). The only remaining time dependence arises from the evolution of the 
electron window function. Its effect is to eliminate the contribution from scales larger than several 
Tg. Since Vg decreases with redshift, its evolution moves the peak of the bias to smaller scales, which 
should have higher bias. The bias approaches a constant (~ 6 for all three cosmological models) 
at large scale, but it drops quickly at very small scale (due to the gas window function). This 
result is consistent with Refregier et al. (1999). One point worthy of mention is that, the behavior 
that bp approach a constant at large scale does not validate the constant bias model, which always 
produces a unit pressure-density cross correlation coefficient. But as we will see in Sec. 4, this 
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is not the case. The pressure fluctuation variance peaks roughly at the scale of cluster core radii 
(~ 0.13h~^ Mpc), which is the direct result of the gas window function (r^ ~ 1/Sh~^ Mpc). 



3. Statistics of Sunyaev-Zel'dovich effect on CMB anisotropies 

The SZ effect on CMB anisotropies is fully described by the n-point correlation functions 
between temperature distortion 6 in different sky directions. The CMB anisotropy is decomposed 
into the harmonic components by: 

e(g) = ^(g) = f2 o^imYUq) (11) 

The CMB power spectrum is defined as Ci = aimaim)/^^^ + 1)- The bispectrum is defined 

as B^fl^^"^'^ = (a™^ a[^^ aj^^ ) . They have been discussed by many authors as mentioned in the 
introduction section. Our model successfully reproduces these results. In contrast to the Press- 
Schechter models, our model relies only on direct statistical correlators that have been proposed 
semi-analytically with only a few free dimensionless parameter, the Q saturation values, that have 
been verified in simulations. In the Press-Schechter model, the shape of the correlation function 
depends directly on the assumed structure of halos, which is a free function of the model and 
empirically measured in simulations. The radially averaged structure of halos, however, is not 
necessarily a direct predictor of correlations, unless one assumes halos to have no substructure. It 
is reassuring to see that our approach, which originates from the diametrically opposite theoretical 
bases from Press-Schechter formalism, produces consistent results and agrees with simulations. 



3.1. SZ power spectrum 

Since we are interested in small angular scale where the line-of-sight CMB survey depth is 
much larger than any correlation length, we adopt Limber's equation (Peacock 1999). In the 
Rayleigh-Jeans region we get the SZ power spectrum: 

^' = (2Z + 1)3 I ^sz{k, z)\k=i/,C{x)xiz)dx{z) (12) 

where 

The results are shown in figure 4 and 5. We find that (1) The SZ effect exceeds the primary CMB 
anisotropy at I ~ 2000. (2) The Q of SCDM is much smaller than that of OCDM and ACDM. 
This is due to the smaller as and faster drop of the gas temperature in SCDM with increasing 
redshift. (3) Because P{k) oc cr| ^ (2 in linear regime and 3 in stable clustering regime) and 
Pp{k) (X B4 oc P^{k), we at once find that Ci oc Ug ^. (4) The peak contribution to C; depends 
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on the angular scale, varying from z ~ 0.4 at / ~ 1000 Xo z 1.5 for I ~ 10000. Smaller angles 
probe the more distant universe. These results are consistent with existing simulations and Press- 
Schechter formalism results. 



3.2. SZ bispectrum 

The SZ bispectrum is the projection of the pressure bispectrum B^{ki,k2, fca) defined by 

(yp(ki)2/p(k2)2/p(k3)) = Sf(ki,k2,k3)(5ij(ki +k2 +k3) (14) 

Since oc 5^, is related to 6-point density correlation function and in principle can be computed 
using our model. But due to the complexity of the 6-point hierarchy and the absence of extensive 
numerical tests, we adopt the approach of Cooray et al. (2000a) and assume a simple pressure bias 
model where 6p{k,z) = bp{k, z)6k{z) with bp{k,z) = Pp{k, z)/ P{k, z). The pressure bias model 
states that the pressure is a linear convolution over the density. Wc note that this model builds 
in certain assumptions, which are invalid for treating cross-correlation issues address in section 4. 
It presumably provides a quick order of magnitude estimator, but fails on qualitative issues like 
galaxy-SZ cross-correlation coefficients. 

In the bias model the pressure bispectrum is given as 

3 



Sf(ki,k2,k3;z) = ^^(/ci, z)53(ki, k2, k3; z). 
The SZ skewness parameter, 



63 = 



((G-e) 



(15) 



(16) 



((0-6)2)2 

is the easiest observable projection of the SZ bispectrum, so we show its computation as an example. 

Instead of working in multipole space, we derive following equations adopting the same ap- 
proximation as Limber's equation: 



27r 



(TTPeiz) 

nieC^il + z) 



C{x)dx{z) X 



bl{k,z)P{k,z)kdk 



(17) 



(0(n)3) 



0^ 



i "'.'/.^'\0 'C{x)dxiz) 



(27r)2 J ^mec^{l + z)' 

3 

bp{ki, 2;)53(ki, k2, k3; z)kidkik2dk2 



i=l 



(18) 



where k3 = ki + k2 and ki || k2. 6 = —2y ~ 5 x 10~^ is the mean temperature SZ distortion. 6^ 
and 0^ are from the corresponding Gaussian term in the correlation function, as explained in Eq. 
(8). As a reminder, C{x) is the geometric function defined in Eq. (3). Though the unsmoothed 63 
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is not directly observable due to the absence of the window function, from the theory viewpoint, 
it represents the raw non-gaussianity properties. Our calculations find that 63 ~ 10^ (Table 1). 
We defined the skewncss parameter (16) as the ratio of differential measurements, which might 
result from interferometers in the Raleigh-Jeans regime. For comparison, Cooray et al. (2000a) and 
Cooray (2000b) considered the absolute temperature distortion and defined the skewness parameter 
as With multi-frequency information the absolute y could be measured at each point 

in space, and that definition of skewness could in principle be observed. The results are not too 
different: In our model, (0^)/(G^)^ ~ 10^, which is similar to the result of (Cooray 2000b) with 
the maximum virialized mass ~ IO^^Mq. Further work considering window function filtering and 
the general calculation of the SZ bispectrum can be done following the method of Cooray et al. 



In the SZ effect, the redshift distribution of intergalactic gas is lost. Taking advantage of the 
cross correlations with other surveys having redshift information, we may be able to statistically 
extract IGM 3-D correlations and evolution. Because the SZ effect is mostly contributed by non- 
linear structures at z ~ 0.5 — 2, it should have a strong cross correlation with galaxies at that 
redshift range. The general idea is to use a galaxy survey with coarse redshift information, where 
the distance information only needs enough accuracy to resolve the time evolution of the correlation 
function, for example ~ 0.2, which can be reached by photometric redshift surveys. One then 
correlates each redshift bin with the SZ map, and obtains the relative contribution of that SZ slice 
to the total projected map. 

Several questions that pose themselves in this procedure are: How big a survey area and depth 
does one need? What fraction of the total SZ fluctuations can be accounted using cross-correlations? 
Are there optimal procedures of cross-correlating? Since SZ fluctuations are a non-linear function 
of the galaxy fleld, what fraction of the signal will ever be accounted for by this approach? Do we 
expect the SZ-galaxy cross-correlation coefficient to depend on redshift? How big an uncertainty 
might arise from a time and scale dependent galaxy bias? In this section we will address each of 
these questions using our model, and obtain quantitative estimates. 

We assume a linear bias model for the galaxy number overdensity 5g = bS. The bias b is taken 
to be a constant for simplicity. Then, the projected galaxies number overdensity can be expressed 



4>{z) is the selection function (which we will want to vary a posteriori using the photometric redshift 
information), zg is the redshift survey depth and C(x) is the geometric function defined before. 



(2000a). 



4. Extracting redshift information from SZ-galaxy cross correlation 



as 




(19) 
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The corresponding multipole moments are: 

= (2[Tl)3 io ^)\k=l/x4>\z)C{z)x{z)dx{z) (20) 

where = 5^A^. The multipole moments of SZ-galaxy cross correlation are: 

Cf^'"" = ^^'''^ ^lz,Gik, z)\k=i/A{z)C{z)x{z)dxiz) (21) 

with A|^ (j(fc,2;) = ^^e^li^^) 2:^fe^^p,G(fc)- ^'p.G is the Fourier transform of (yp(x)(5G(x + r)). 

Pik)fe{k) + 7^1 Bsih, k2, k)Wg{h)fe{k2)d^k2 

P„G(k,z) = bx (22) 

We have used k = ki + k2 . It is useful to define the cross correlation coefficient 

, ^ ^f^«<f-'> ^ , , ^'■■<'<*-^) (23) 

In -Pp,G; the i?3 term is dominant. In ^4 there are 16 hierarchical terms and in there are three 
hierarchical terms. Different term dominates in different regions. Calculation shows that, roughly 
there are three regions: (a) k < O.lh/Mpc. 4 terms ( Pi2(2PiP4 + -P1-P3 + P2P4) ) dominate B4 
and 2 terms ( P3P1 + P3P2 ) dominate B3. r ~ Q^/Q]/'^ 0.9 (n 1.5)^. (b) A; < O.lh/Mpc < 

lOh/Mpc. Each hierarchical term in i?3 and B4 has about the same contribution to -fp,G and Pp, 

1/2 

respectively. r(/c,z) ~ 3/4 x Q3/Q4' ~ 0.7 (n ~ —1.5). This region contributes most of the SZ 
effect, as seen from figure 3. Unless explicitly notified, hereafter we will adopt the value of r in this 
region. We only show the result of this region in Figure 6. (c) k > lOh/Mpc. This is the opposite 
case to the case (a), r ~ Qs/V^'^Qa — 0.3 (n ~ —1.5). No significant time dependence is found. 
The corresponding cross correlation coefficients in multipole space are: 

Corr{l, (f)) = 24 

[Cszii)CGiiA)r 

Wc now use the photometric redshift information from the galaxy survey to vary the selection 
function 0. We pick the weighting which maximizes (24). This allows the measurement of one 
number in the SZ to yield a full function of redshift z. Since (24) depends on one variable /, we can 
in principle measure an optimal redshift weighting (piz) at each I. The cross correlation variation 
has allowed us to measure a two dimensional cross-correlation function from a one dimensional 
observable in the SZ and a two dimensional observable in the galaxies. 

^Though k is small, contribution to Pp{k, z) and Pp^aik, z) are mostly from nonlinear regimes with — 2 < n < — 1. 
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The selection function to maximize Corr{l) is obtained from the variation ^'^"^^^''t') = q. We 
denote this selection function as <Pm{1,z) and find that 



^m{1,z) = a{l) 



(25) 



Here, a{l) is a constant to be determined later by the observational data. The optimized C^^orr 



is: 



C^orr{l) 



It"'" ^lzC{x)x{dx/dz)dz 



1 1/2 



/o" ^C{x)x{dxldz)dz 
jzcMB A|^C(x)x(da;/dz)dz 



1/2 



(26) 



Results are shown in figure 7. Following the same estimation as in r{k,z), C orril) c± 3/4 x 
Qz/Q^J^ ~ 0.7. The observed C^orr{l) may be smaller because of the limited galaxy survey 
depth. 

The observationalroccdurc to extract the rcdshift information is as follows. 1. Start with a 
random guess for (t){z), e.g. (p = 1 . 2. Given a photometric galaxy survey and SZ survey, measure 
Csz{l), Cg{1,4') and Csz,g{1,4') ^om angular correlation functions {Q{h)@{h + 9)), {^{n)^{h + 
§)) and (0(n)^(n + §)), respectively. 3. Vary cf) to maximize Corr{l) for a specific Z, therefore 
obtain (j)M{l,z). 4. Apply Eq.(25) to infer ^(//^(z), z), up to a constant l/a{l), from the 
directly measured from the galaxy survey. Furthermore, Eq. (21) enables us to infer a{l) from 
the observed Cf^^f . 



ail) 



2 X IGtt^ 



4'm{Ij z)Aq[1/x{z), z\C {z)x{z)dx / dzdz 



I '^l,obs 



(27) 



5. Combine Eq.(23) and Eq.(25) to obtain the time resolved SZ power spectrum, up to a factor 
r{l/x{z),zr^. 

Alz{k,z) = ct>l,{l,z)Al{l/x{z),z)/[r\l/x,z)a\l)] (28) 

In our model r{l/x{z),z) is almost a constant 2± 1/0.7^ 2± 2 independent of cosmological models 

or redshifts (sec Figure 6). This gives a theoretical estimate of A'g^{l/x, z). 6. nomalizes 
the observed SZ anisotropy multipolc C°^* and the SZ temperature variance. Apply Eq. (12) to 
integrate A'g^ obtained above with 1/r^ = 1 and compare with (7;°''^, we will obtain the averaged 



r'^{l/x{z),z) ' 



-^r^^ J ct>M{l,z)Al[l/x{z),z\/a\l)C{x)x{z)dx/dzdz (29) 



7. Follow the same steps to obtain A^g^ g{^/x{z),z) and A^g2{^/x{z),z) at different I. 

Since for Csz{l) at different angular scale I, {l/r'^{l/x, z)) is determined roughly by r{l/x{zp), Zp) 
[zp ~ 1 is the redshift with peak contribution to Q. See Fig. 5), we can even get some idea about 
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the scale dependence of r{k,z). Thus, the total projected SZ autocorrelation give a consistency 
check on the reconstructed time resolved power spectrum from the galaxy-SZ cross correlation. 
Furthermore, the galaxy bias and its time dependence have completely dropped out of the calcula- 
tion, and are thus not expected to affect the results at all. Then, in principle, SZ-galaxy correlation 
plus SZ CMB anisotropy provide a consistent and powerful method to extract all time evolution 
information of the IGM pressure power spectrum. 

Noise and cosmic variance put constraints on the feasibility of our procedure. (1) Limita- 
tion of CMB resolution degrades our method. The measured range of k is [li/x{z),l2/x{z)] ~ 
[li/3000z, Z2/30002;]/i/Mpc. Here, [hjh] is the range of the CMB experiment. In order to detect 
the peak of ( around k = 3/i/Mpc as shown in Fig. 3), z < h/QOOO. For CBI (630 < I < 3500), 
we are only able to detect z < 0.4. AMIBA will measure I < 28500 and South Pole Submillimeter 
Telescope (2003) will measure / < 40000. They will allow us to measure the gas power spectrum 
up to 2; ~ 3 and z ~ 4, respectively. (2) Observational errors impose further constraints. Suppose 
that the galaxy survey covers a fraction fa of the sky and the i-th survey region (For example, if 
the redshift accuracy is Az, then we can divide the galaxies into redshift bins with < z < Az, 
Az < z < 2Az, etc. ) have Ni observed galaxies. The CMB observation covers a fraction fcmb 
of sky and the Ci is averaged over the band [/ — Al/2,l + A//2]. Then the galaxy number count 
causes the Poisson error: 

^ . ^sz^ ^ [fa X rmnimr^/^ (30) 

(3) The cosmic variance of the C; also cause errors. Recalling that Ci = J2 '^imO-im/ (^^ + -*-) 
aim oa 6T (X bp6, we get: 



ACi_ UCf)-{Ci)^ S,al 



Ci V {CiY V (2^ + l)A//cm6 V 10(V2000)AZ/e^6 



(31) 



Here, a\ is the density dispersion over smoothing scale R ~ /i/Mpc. We already use the typical 
value of 54 ~ 40 and cf\{z = 1) ~ 10. The corresponding error caused in 0m is: 




2 [(10(Z/2000) A//^,)-i + X minm-'] (32) 

M 5<l>'^ ''I'M 

Recalling that A|^ oc requiring a 40% accuracy on A|^ would impose that (a) fc x min{Ni) > 
10^. Each survey regions must be large enough in order to contain sufficient number of galaxies 
and must be small enough to ignore the evolution. We may choose redshift bands of each survey 
region Azi ~ 0.1. Then the number of galaxies observed Nq has to satisfy Nq > W^zc/Azi/fG- 
For SDSS (Sloan Digital Sky Survey (2000)), which covers one quarter of the sky and probes more 
than one million galaxies with photometric redshift up to 2; ~ 1, the requirement Nq > 10^ is easily 
satisfied. The measurement of the intergalactic gas at redshift z requires the galaxy survey at least 
up to that redshift (equation 25), so we need deeper galaxy survey in order to probe the gas beyond 
z ~ 1. (b) fcmb^l > 10^(2000//). For CMB experiments with relatively lower resolution, larger 
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sky coverage is required. For example, though Planck only measures I < 2000, it covers the whole 
sky and therefore satisfies this condition. For those with much higher resolution such as AMIBA 
and Submillimeter Telescope, the required sky coverage can be relaxed to the order of 1%. 

This variation method only depends on the assumption that the cross correlation coefficient 
r(k, z) is approximately a constant, which is the direct result of the hierarchical model and has only 
weak dependence on the gas model and cosmologies. Furthermore, the hierarchical model is strongly 
supported by the consistency of the CMB SZ power spectrum dependence on erg and the behavior 
of the gas bias between our model and simulations. Since the averaged r{k, z) is measurable, our 
method does not rely much on the theoretical value of r(A;, z) and thus observationally consistent. 



5. Non-gravitational heating effect 

Pen (1999) has shown that the IGM has most likely been preheated by non-gravitational energy 
sources with energy injection Ejqc ^ 1 keV per nucleon. This section is devoted to consider this 
effect. Because the relation between the gravitational heating and the non-gravitational heating 
is very uncertain, we only consider two extreme cases. The first one is that the non-gravitational 
heating is perfectly correlated with the gravitational heating, then we can change Eq. (5) to: 

fcT, = i(l + ^)^^(*-t) (33) 

We will append all former results from gravitational heating with a superscript 'A' (Adiabatic). 
Here, 13 = Tj^c/Tg^ represents the ratio of the non-gravitational heating and the gravitational 
heating. kTj^o = 8/3/(3 -|- 5X)Eng- Enq ~ 1 keV corresponds to /? ~ 1. All former results are 

_ _ _ 2 — 

not affected by this change except Tg, y oc Tg, Ci oc Tg and 63 oc 1/Tg due to the dependence 
Tg (X {1 + (3). Figure 1 and 4 need to be changed correspondingly. Figure 5 and 7 will change only 
when (3 is time-dependent. All other figures remain the same. 

The second case is that non-gravitational heating is uncorrelated with local density, then, Eq. 
(5) changes to: 

= ^sT^^* " + ^^^^^^^ ^^^^ 

Tg = Tg^+TNG = [^+P{z)]Tg'^- Then, the corresponding new results following the same definitions 



are: 



Pp{k, z) = [P/(fe, z) + (27r)3fe(k)(2/? + /32) + 2pPiz,5{k, z) + P^P{k, z)] /(I + (3f 

= l^'f ^ X Pik, z) (when k ^ 0) (35) 

Pv,G{k, ^) = X [P^^,{k, z) + (3P{K z)] = 6 X !:l^|±^P(fc, z) (36) 
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bp{k, z) = (37) 

rik, z) = . '''^^ (38) 



Here, P^^ = P^qI^ is the Fourier transform of (?/^(x)(5(x + r)). The expression of the optimal 
selection function <\)m (Eq. 25) is not affected at all. Equation (38) tells us that < r <\. As 
expected, when /? ^ 6^ , we reduce to our former results and when /3 » 6^ , we obtain the bias 
model. There are three regions: (a) k S> l/r^ ~ 3/t/Mpc. We always have /3 S> 6^ , — /3/(l + /3) 
and r — > 1. Our redshift extraction method works well with r ~ 1 in this region. Because this region 
is dominated by non-gravitational heating, high resolution CMB experiments are able to measure 
A^^ in this region and enable us to directly obtain the information of the non-gravitational heating. 
Since 6p does not approach zero as before, A|.^ oc increases with increasing k. This behavior 
will move the peak of SZ power spectrum to larger / while increasing the amplitude of SZ power 
spectrum, (b) A; < 6/i/Mpc. Then, 6^ ~ 8 S> /3(0) ~ 1, so hpiz) will decrease roughly by a 
factor (1 -|- (3{z)). The resulting r is bigger than our former results, but it still remains roughly a 
constant with respective to space and time, so our redshift extraction method works well, (c) The 
intermediate region ( the range ~ /?, roughly 6h/Mpc ^ k < lOh/Mpc), r varies from ~ 0.7 
to 1. Our redshift deprojection method still works, but with a larger error. The discussion of the 
observational requirement is not affected at all. The new Corr should differ only a little bit from 
the former result due to its dimentionless definition. Roughly, C; oc (1-1-/3)^, bispectrum oc {l + (3)^ 
and skewness @s oc {1 + . 

Our knowledge on the non-gravitational heating is very limited, and presumably depends on 
the poorly understood physics of star formation and supernovae dynamics. It is not even known if 
the heating was pre or post structure formation. We have surveyed two simple but different models 
of non-gravitational heating to demonstrate the range of effects it may have on the correlations. The 
above equations show the basic procedure to include the non-gravitational heating and the possible 
effects. So we do not plan to constrain to a highly hypothetical model and do the calculation. For 
the purpose of estimation, we may either assume a step function for Tjvg {Tng{z) = Tj^ai^) for 
z < zng- Otherwise, Tng = 0.) or Tng{z) = T]vg(0)(1 + z)~^ . The upper limit of 'y' parameter 
y < 1.5 X 10~^ from COBE/FIRAS (Fixsen, et al. 1996) put constraints on the value of zjvG and 
7. For example, when T/vc = 1 kcV, 7 > 1.5 or zj^g < 3. Further simplification can be made by 
taking the fitting formula of {bp)^ and ~ f"^ ~ 0.7. We show some examples (Figure 8,9 and 
10). 

When combining CMB experiments with galaxy surveys, we have four observables: Q, Ag(fc, z), 
^"szoik-i ^) = (f'Mik, 2;)Ag(fc, z) and f{k). If r(k, z) is not a strong function of time and scale, we 
can solve for /\^g^{k, z). If the behavior of bp{k, z) as discussed above is correct, we obtain the fifth 
observable (3{z) from the small scale behavior of Al^- 



6. Discussion 



Let us now address some potential systematic shortcomings in our simplified model. (1) Our 
model relies on the adopted gas-dark matter density relation (Eq. 4) , the (non-local) temperature- 
density relation (Eq. 5) and the non-gravitational heating. In case that only the gravitational 
heating is included, apart from intrinsic problems addressed in section 2, our model contains two 
free parameters: the gas-dark matter smoothing length, and the temperature-potential smoothing 
length. These parameters only weakly affect the pressure power spectrum and pressure bias, which 
are defined in terms of dimensionless functions. The normalization is expressed in terms of the 
mass weighted temperature, Tg. We see in Table 1 that the gas smoothing only has a small 
effect. A preheating of 1 keV as proposed by Pen (1999) has a factor of order unity impact on the 
temperature, and about factor of 4 effect on the CMB anisotropy spectrum and 8 on the bispectrum, 
while the skewness will decrease by about 50%. We provide a convenient method to include such 
effects. Hydro simulations arc able to measure the gas-DM relation and the T — ^ relation and 
will further improve our work. Fortunately, the main goal of our model — the extraction of the 3-D 
gas information is least affected, since it depends only on the cross-correlation coefficient (23) and 
Figure 6, which does not vary a lot even including non-gravitational heating. (2) The galaxy bias 
model we adopt may be too simple. In our redshift extraction method of the gas pressure power 
spectrum, the only bias dependency comes from r(k,z). Wc speculate that, even for a realistic 
galaxy bias model, r{k,z) should be close to a constant, which still enable us to extract Al^- (3) 
Ql^*. The value of Q^* in HEPT is a function of the power index of a power law power spectrum. 
We have extended it to the CDM power spectrum and choose the power index at Y^ki/N. Because 
SZ effect is mainly contributed by the non-linear regime with n ~ —1.5, HEPT is applicable and 

1/2 

the resulting Qjv does not vary a lot. More importantly, Q3/Q4 only has a weak dependence on 
n. This behavior ensures that r{k,z) has least dependence on our assumptions about Ql^*- 



7. Conclusions 



We have presented a new tool to compute power spectra and other statistics of SZ fluctuations 
based on hierarchical clustering and scaling. This approach describes the two point correlation of 
non-linear observables, such as the SZ y parameter in terms of two point correlation functions, which 
directly maps to the observed angular power spectrum. We have demonstrated that this approach is 
feasible, and produces results consistent with simulations and Press-Schechter approaches. We then 
addressed the problem of measuring the redshift evolution of the SZ effect through cross correlation 
with galaxy surveys containing coarse photometric redshift information. We have shown that a 
variational method allows the redshift deconvolution which does not depend on galaxy biasing. 
The only model dependent quantity is the cross correlation coefficient r between galaxies and gas 
pressure, whose value has least dependence on the gas density model and the gas temperature 
model. 
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Our quantitative estimates suggest that, when combining CMB experiments either with high 
resohition such as AMIBA and South Pole Submihimeter Telescope or with high sky coverage 
such as Planck and deep broad photometric galaxy surveys such as Sloan, our method is able to 
extract the redshift evolution information of intergalactic gas, such as the full time resolved gas 
pressure power spectrum, even without requiring the knowledge of galaxy bias. It is even capable of 
disentangling the contribution from the gravitational heating and those from the non-gravitational 
heating. This method serves as a powerful probe to this primary component of cosmic baryons 
to high redshift. The model also provides an alternative to simulations and the Press-Schechter 
formalism to calculate the CMB SZ power spectrum and bispectrum. We successfully reproduce, 
when no non-gravitational heating presents, the mass weighted gas temperature Tg{^ 0.35 keV), the 
pressure power spectrum, the pressure bias (~ 8, but scale dependent), the mean SZ temperature 
distortion (~ 5 x 10~^ K), the SZ power spectrum and the skewness parameter (~ 10^). With our 
transform formulas, non-gravitational heating is easily included and we have estimated its effects 
on the correlations. 
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Fig. 1. — The evolution of the gas density averaged temperature. Smaher Tg produces higher Tg. 
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Fig. 2. — Gas pressure bias in ACDM model. Results for OCDM and SCDM are similiar. 
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Fig. 4. — CMB anisotropy caused by SZ effect in various models. Vg = l/3h ^Mpc. Smaller rg will 
move peaks to larger I ( smaller angle) and larger amplitude. 
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contribution to SZ power spectrum from various z 
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Fig. 5.— SZ effect contribution from different redsfiift, defined as r^.f'Jit'o '"^'^"^^^1''?.^?^. . ■ We 

sflow a ACDM model (r^ = ^/h Mpc). Results for SCDM and OCDM are similiar. Smaller r„ 
moves peaks to smaller z and lower height. 
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Fig. 6. — SZ-Galaxy cross correlation coefficients in real space for a ACDM model. Results have 
almost no dependence on cosmological models and gas model parameters. The region shown is the 
most relevant region of SZ effect, as shown in figure 3. The behavior of r{k, z) in other regions is 
estimated in Sec. 4. 
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Fig. 7. — SZ-Galaxy cross correlation coefficients in multipole space. Results have almost no 
dependence on cosmological models and gas model parameters. 
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Fig. 8. — bp(k,z) and Ap{k,z) when non-gravitation heating is considered. kT^c = 0.39keV(l + 
z)~^ is assumed. The time evolution of bp is partly due to the evolution of (3. 
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Fig. 9. — Ap{k,z) when non-gravitation heating is considered. kT^c = 0.39keV(l + z)^^ is as- 
sumed. The distinctive behavior of at small scales {k > lOh/Mpc) is due to the non-gravitational 
heating and the measurement of this region will help us to extract the non-gravitational heating 
information. 
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Fig. 10.— 
8. 



r{k,z) when non-gravitational heating is considered with same assumptions as in figure 



